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5 . METHOD FOR MONITORING WOOD-DRYING KILN 

STATE 

TeclMiical Field 

The invention relates to control systems for kilns used to dry wood. 

10 

Background 

Wood naturally contains moisture. However, for wood lumber to be useful for 
building, the moisture content must fall within certain constraints. Too much or too little 
water in the liunber can cause serious problems in the subsequent processing and use of 
1 5 the lumber, including warping or difficulties with painting or gluing. In order to provide 
a useful product, most lumber mills dry their lumber under controlled processes, such as 
in a kihi, prior to sale. It is difficult, however, to ensure that every board in a given 
charge will have exactly the same moisture content after drying. 

Part of the problem Ues in the nature of the wood itself; wood is a non- 
20 homogenous material that does not dry at a constant rate. Even within a charge with all 
boards of the same size and same type of wood there are many factors that can affect 
drying rate. Individual different boards may have natural variability in their drying rate, 
may start out with different moisture content, may be composed of different percentages 
of sapwood and hardwood, or may contain wet pockets. Moisture variations can also 
25 result from variable drying conditions between different kilns at the same mill or within a 
single kihi charge due to tiie locations of particular boards within a kiln. While the 
variation between charges or kites can be monitored using averages, the variability within 
a charge is best monitored using standard deviations. One method by which this 
monitoring can be done is by using statistical process control (SPC). 
30 SPC works to continuously improve the quality of a process through the constant 

apphcation of statistical methods to data produced at various stages of the process. 
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5 Ideally, disturbances within the process are quickly discerned, then, using the data set, 
causes are assigned to those disturbances. In the best case, this, is done with sufficient 
speed that investigation of the problem, diagnosis, and corrective action can be taken 
before many noncomforming products reach the final state of the process. This is usually 
done through the use of parametric control charts for variables. Generally, by plotting the 

10 average and standard deviation of multiple samples on a control chart and comparing 
those values against computed acceptable limits, the chart can be used to determine if a 
process is going out of control. 

The starting point for SPC analysis is an understanding of the underlying 
distribution of the data under analysis. The conventional assumption has been that wood 

1 5 moisture content data is approximately Normally distributed. However, wood is a non- 
homogenous material, and evm when it is kihi-dried under controlled conditions, the 
moisture content retains large variabiHty through the wood pieces, so a board with an 
acceptable average rating may be quite dry in one section and quite wet in another. 
Additionally, because kiln-dried wood cannot be dried past the point at which it reaches 

20 equilibrium with the air in the Idhi, there is an effective lower limit on moisture content. 
Thus, a Normal distribution has proven less than ideal for SPC analysis of wood drying 
and has made the task of creating a proper model for SPC analysis of wood drying 
di£5cult.It is desirable to design control charts specifically to monitor wood-drying, using 
a better-performing distribution model, if possible. 

25 The theory of variability is the basis for the construction and use of control charts. 

The purpose of any monitoring process is to sort out meaningfiil signals fi-om background 
noise. A certain amount of variability in the moisture content of the dried lumber- 
background noise— exists in the kibi-drying process. Boards dry at different rates due to 
inherent properties in the wood itself, as mentioned above. This natural variability can be 

30 considered a system of "common causes." A drying process that is operating with only 
common causes of variation is considered to be in statistical control. 

A major objective of the control charts is to sort out a meaningful signal firom the 
background noise. When the moisture variability of the dried wood exceeds the 
background noise variability, this might indicate that there are problems associated with 
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5 the kiln itself. This variability in moisture content usually arises from three sources: 
improperly adjusted kibi equipment, operator errors, or defective raw material. Such 
variability may be large when compared to the background noise, and it usually 
represents an unacceptable level of process performance. These sources of variability are 
called "assignable causes." A process that is operating in the presence of such causes is 

10 said to be "out of control." A kiki can operate in the "in control" state for a long time, 
producing acceptable moisture content levels in lumber. Eventually, however, assignable 
causes will occur, resulting in a shift to an out of control state, where the lumber 
produced no longer conforms to requirements. 

A major objective of the control charts is to quickly detect the occurrences, or the 

15 trends toward assignable causes of process shifts, which allows investigation mto causes, 
and corrective action to be undertaken before many nonconforming units are produced. 
Control charts can also be used to estimate the parameters of the kibi-drying process, to 
determine process capability, and to incrementally improve the process. In general, 
control charts are an effective tool that can be used to reduce process variability. 

20 Many quality control methods suggested in the literature are concerned with the 

average moisture content alone, and therefore use methods such as sample estimation, or 
go/no go decision criteria, such as acceptance sampling. Other methods have been tried 
that do concem themselves with variation, but are based on the assumption that moisture 
content has a Normal probability distribution. This assumption tends to erroneously 

25 consider certain in-control conditions as defects or out-of-control conditions. 

Furthermore, assumptions of Normality underlie many other common statistics, such as 
correlation coefficients, t-tests, and others that may be used in standard statistical 
analysis. The failure to understand that a distribution is not Normal leads to a rippling 
statistical effect, with invalid assumptions producing incorrect answers which are 

30 themselves imperfectly analyzed, and so on. 

Lognormal, rather than Normal, probabihty plotting for kihi analysis has been 
proposed by "Moisture Control During Kihi Drying," McMahon, 1961, in conjunction 
with frequency distribution analysis and quality control methods. These techniques were 
used to estimate the average moisture content and the percentage of moisture content 
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5 outside of certain boundaries. McMahon's methods were only 25)proximate, however and 
introduced procedural errors that. were potentially large. Additionally, McMahon's 
technique failed to properly account for the lower moisture content limit of kihi-dried 
wood. What is needed is a method that accounts for the lower moisture content limit that 
kiln-dried wood approaches and monitors kiln conditions more closely. 

10 

Smimmary 

A method for monitoring the drying of wood in a kiln is described herein. The 
method uses a set of known good samples of moisture content data for wood dried in a 
kiln to estimate a Lognormal statistical distribution which approximates the known good 

15 set of data. The method estimates three parameters of a Lognormal distribution: the 
threshold parameter, the shape parameter, and the scale parameter. Then one or more 
new sample of moisture content data is obtamed from a subsequent drying in the kihi. 
The method develops one or more control charts containing center lines and control 
limits which monitor at least one of the estimated parameters. Next, data from the new 

20 sample or samples of moisture content data is transformed to represent the parameters on 
the control charts. Finally, a determination is made about whether the kihi is operating in 
an in-control or out-of-control state. 

A computer-readable medium is also described herein. The computer-readable 
medium contains instructions which, when executed, obtain a known good set of samples 

25 of moisture content data from wood dried in a kihi. Then the instructions use the known 
good data set to estimate parameters for a Lognormal distribution which approximates the 
known good data. Next the instructions obtain new moisture data for wood dried in the 
kiln. The instructions then create control charts based on the estimated Lognormal 
parameters, transform the new data to be represented on the control charts, and determine 

30 whether the kiln is operating in an out-of-control state. 

Additional features and advantages of the invention will be made apparent from . 
the following detailed description of embodiments that proceeds with reference to the 
accompanying drawings. 
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5 

Brief Description Of The Drawings 

Figure 1 is a diagram of a suitable environment in which described embodiments 
may be implemented. 

Figure 2 is a diagram of a suitable environment for measuring moisture content of 

10 wood. 

Figure 3 is a diagram showing examples of wood moisture content distribution 
patterns. 

Figure 4 is a diagram showing examples of different distribution models. 
Figure 5 is a block diagram of a method to develop control charts 
15 Figure 6 is a block diagram of a method to estimate Lognormal shape and scale 

parameters. 

Figure 7 is a block diagram of a method to create a control chart. 

Figure 8 is a block diagram of a niethod to sample, transform and plot sample data 
to determine if a kiln is out of control. 
20 Figure 9 is a diagram showing examples of control charts produced from a 

Lognormal distribution model. 

Detailed Description 

The following description is directed to improvements in creating quality control 
25 charts so that they more accurately and realistically monitor the variability of moisture 
content in kihi-dried lumber. These control charts are then used with more-recently 
collected data to determine if a kihi is out-of-control. 

1. Moisture Content Measuring Environments 

In an exemplary embodiment, moisture control measurements are collected from 
30 kiln-dried lumber. Figure 1 depicts one environment for measuring moisture content in 
wood. In the illustrated embodiment, a lineal moisture content sensor 102 is controlled 
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5 by a moisture content control console 104. By adjusting various settings on the console 
104, the sensor 102 can be made to measure the moisture content of varioxis types, sizes, 
and shapes of wood. When data is received by the console 104, the data is sent to a 
computer 106, nmning moisture content monitoring software 108 (not pictured). As will 
be discussed below, the data received by the computer 106 can be analyzed by software 

10 1 08 to determine if any out-of-control situations exist in the drying process. 

With respect to Figure 2, an environment featuring various moisture sensors in a 
lumber mill is shown. In the illustrated embodiment, wood boards 201, which have been 
dried in a kihi, are sent through various in-line continuous moisture sensors, including a 
lineal moisture sensor 204 and a transverse moisture sensor 206, both controlled by 

1 5 moisture control consoles 202. Also illustrated is computer 208, running moisture 
content analysis software 210, which collects data from the moisture content control 
consoles 202 for analysis. In another embodiment, moisture content readings may be 
made by types of sensors other than those illustrated, including a hand-held sensor, a pin- 
type sensor, or a contact sensor. In yet another embodiment the computer 208 and 

20 software 210 may be integrated into individual control consoles 202 or the system may 
incorporate a single master computer which communicates with every control console 
and sensor to receive and analyze moisture content data. 

2. Particulars of Wood Moisture Content Analysis and Lognormal 
Distribution 

25 As was mentioned above, the unique nature of wood makes it difficult to dry 

every board in a kiln charge to the same moisture content, providing a variation. Even 
when boards are made of the same type of wood, are cut to the same specifications, and 
are dried for identical times and at identical temperatures, variations can still occur 
between kibis, between charges within a kibi, between boards in a single charge, and 

30 even within a single board. Thus, in order to provide as accurate a measure as possible of 
the stability and consistency of the drying process, statistical methods can be used. ^ 
However, because the data concerns moisture content in wood, additional consideration 
must be given to the statistical model used to analyze the data. 
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5 Figures 3a-d illustrate various examples of distribution for various samples of 

wood moisture content data. The examples shown were produced by the inventors icom 
studies done by Pratt in 1953, McMahon in 1961, Bassett in 1969, and Milota and Wu in 
1995. For ease of description, relevant aspects of the moisture content data distribution 
examples will be discussed with reference to the histogram displayed in Figure 3d. 

10 Of chief concern for statistical modeling are the finite left end and positive 

skewness of the distribution. As the example histogram of Figure 3d demonstrates, and 
as the histograms of the other figures support, there is a tendency for moisture content 
distribution to drop oflF dramatically on the left end of the distribution. This is illustrated 
in Figure 3d by the location 304. Typically, this dramatic drop-off is due to the fact that 

1 5 the wood cannot have a moisture content lower than the eqmlibrium moisture content 
(EMC) set up in the air of the kiln during the drying process, and in any case can never 
have a moisture content lower than zero. Thus there is a physical lower-bound to the 
moisture content that can be achieved by kihi-drying, and so a lower-bound to any 
moisture-content data that can be measured. Additionally, because not every piece of 

20 wood dries at the same rate, moisture content data fi-om kihi-dried wood exhibits a . 

positive skewness, as demonstrated by the region 306. This means moisture content data 
distributions tend to show a trailing right-hand side even while the left had side has a 
definite edge. While the illustrated characteristics of moisture content data distribution 
are not necessarily found in every sample of dried wood, they are common to most 

25 samples taken, and are thus important to consider when searching for a model that fits 
well. The system and methods described herein account for these features through the 
use of a 3-parameter Lognormal distribution model. 

Figures 4a and 4b illustrate examples of Normal and Lognormal distribution 
models. In the example graph of Figure 4a a Normal distribution curve is demonstrated. 

30 As the graph of Figure 4a demonstrates, Nomial distribution curves are horizontally 
symmetrical and are positively and negatively imbounded. Thus, the Normal curve is 
fi-equently not the best fit for the types of distributions mentioned abovie with respect to 
Figures 3a-d. 
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5 Lognormal distribution models, however, do exhibit some of the characteristics 

described above. A Lognormal distribution naodel describes a variable that is distributed 
so that the natural logarithm of the variable has a Normal distribution. The graph of 
Figure 4b illustrates a sample Lognormal distribution curve. Unlike the Normal 
distribution curve in Figure 4a, the Lognormal curve demonstrates a positive skewness 

10 similar to that described in region 306 of Figure 3d. The Lognormal distribution curve of 
Figure 4b also demonstrates a left-hand drop-off similar to that illustrated in location 304 
of Figure 3d. In fact, as will be described below, the Lognormal distribution utihzed in 
the present system includes a threshold parameter which provides a bounded minimum to 
the distribution, thus more closely approximating the behavior of kihi-dried wood. 

1 5 The Lognormal distribution model utilized in the present system and method 

makes use of three parameters, a threshold parameter 0, a scale parameter and a shape 
parameter cr. As described above, the use of the threshold parameter is to provide a lower 
bound for the distribution model, i.e. a minimum amount of water that could exist under 
the drying conditions being studied. The scale parameter essentially serves to stretch or 

20 compress the curve on the horizontal axis, while the value of the shape parameter subtly 
affects the shape of the curve. The use of parameters such as these is well understood in 
general statistical distribution analysis, although not with regard to analysis of wood 
moisture content. Mathematically, the general form of this Lognormal probability 
distribution function with three parameters is giyen by the function: 

While data that follows a Lognormal distribution can be described with the 
function (1), the moisture content data is not analyzed directly in the system and methods 
described herein. Traditionally, a control chart plots the arithmetic mean and standard 
deviation (or variance) of multiple samples against a central tendency line and control 
30 . limits in order to monitor process conditions. However, for variables that follow a 
Lognormal distribution, the distribution of standard deviations is not well understood. 
Additionally, because of the nature of Lognormal distributed variables, the arithmetic 
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5 mean of a sample is not as good an indicator of central tendency as the geometric mean. 
Thus, monitoring of averages and. standard deviations of the direct sample data is not as 
useful as other methods. 

Instead, the present methods and system perform analysis on logarithmically- 
transformed data in order to monitor the shape and scale parameters themselves. Because 

10 the logarithm of a Lognormal variable is itself Normally distributed, it is known that for a 
variable which is distributed Lognormal with a threshold parameter 0, a scale 
parameter and a shape parameter <t, then the logarithmically-transformed 
variable Y = hi{X - ^) has a normal distribution with mean /i and standard deviation a. 
Because methods of SPC are well known for Normal variables, techniques for monitoring 

15 the mean and standard deviation of the logarithmically-transformed variable will provide 
information about changes in the shape and scale parameters of the Lognormal 
distribution. It is the monitoring of these parameters that the techniques below focus on. 

3. Creatiora of Control Charts for Wood Moisture Content Data 
Figure 5 illustrates a block diagram which describes, in one embodiment, a 

20 process for creating a control chart for wood moisture content data. Starting at block 
502, known good data is collected so that proper limits and a center line for the control 
charts can be developed from the data. Because differences in wood type or kiln can 
affect the drying process, in the preferred embodiment the collection of known good data 
is done for the specific kiln and under the same conditions in which subsequent drying 

25 will take place. In one embodiment this is done by taking multiple samples, with each 
sample containing multiple wood moisture content readings. These conditions can 
include, in various embodiments: the species of wood, the type of lumber made of the 
wood, the method of stacking wood in the kiln, and the drying schedule. In one 
embodiment the data is divided up for , analysis into subgroups, allowing greater 

30 discernment of kihi drying problems. Thus, analysis can be performed at once for all 

boards in a kiln charge, between packs of boards in a kiln charge, or even between boards 
in a particular pack. This provides flexibility to more readily identify and correct 
problems. 
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5 In one embodiment, in order that the data can be trusted to properly represent 

moisture content data distribution, die data collection is done at a time when the kihi 
being used is known to be "in-control," i.e. operating within normal parameters and 
producing conforming product. In the preferred embodiment, care is taken during the 
process of block 502 to ensure that the known good data taken is done randomly and 

10 independently. In one embodiment this is done through a random sampling with 
replacement, and by ensuring that no two measurements had any overlap. 

Next, at block 504, it is determined if the known good data fits a Lognormal 
distribution. In the preferred embodiment the process of block 504 is performed so that 
in case, for some particular set of conditions, the data does not follow a Logn9rmal 

1 5 distribution, subsequent analysis can follow Normal (or other) distribution techniques 
rather than the Lognormal techniques described herein. In various embodiments, this 
determination is performed through the use of graphical and/or numerical methods, such 
as investigation of symmetry and probability plotting of the Lognormal variable Xand 
goodness-of-fit tests of the Normal distribution to the transformed variable Y, 

20 If, at block 504 it is determined that the known good data follows a Lognormal 

distribution, then at block 508 a Lognormal distribution is estimated that approximates 
the known good data. In the preferred embodiment, the threshold parameter ^is first 
approximated &om the sample at block 510, and then approximate values for the scale 
parameter // and the shape parameter crcan then be calculated fi-om 0 at blocks 512 and 

25 514, respectively. 

With respect to the process of block 5 10, the threshold parameter in the illustrated 
embodiment is estimated fi-om the known good data set in order that the other parameters 
can be calculated. While typically the EMC of the kihi during drying is knovm, that 
value is not substituted for the threshold in the preferred embodiment. This is because, in 

30 order to quickly pull water out of wood, the wood is never allowed to reach the EMC of 
the kihi before the EMC is changed. Thus, in the illustrated embodiment, it is not 
^ assumed that the EMC is the threshold of the distribution. 
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5 In one embodiment, this is done through the technique of using local maximum 

likelihood estimators for as will be understood. In pne embodiment, a local maximum 
likelihood estimator for 9, using an ascending ordered sample x^,x^ ,K , is given by: 

fn \2" 

^ -I (2) 

-"■•?(^)- 

This equation can then be solved using trial and error with linear interpolation to 
10 yield an estimator for the threshold parameter, 6 . In other embodiments, when this 

technique does not yield a vaHd result, an alternative technique using modified moment 
estimators, as demonstrated in ^Modified Moment Estimation for the Three-Parameter 
Lognonnal Distribution," Cohen et al., 1985. Once a sufficient estimate for ^is found, 
estimates for the other two parameters can be found by: 

15 /i = i*2l"(x,-^) (3) 

n J 

d^=l*tln^(x,-^)-[i*|;in(x,-^)J ' (4) 

Relations (2), (3) and (4) come from "Estimating Parameters of Logarithmic- 
Normal Distributions by Maximum Likelihood," A.C. Cohen, 1951. After calculation of 
the Lognormal parameters, the process at block 516 creates one or more control charts. 

20 As was explained above, control charts can measure either the scale parameter, as in 
blocks 518 and 522, or the shape parameter, as in blocks 520 and 522. In the preferred 
embodiment, the shape and scale parameters are monitored, as values for them can be 
determined from the threshold parameter. While creating a control chart that directly 
monitors threshold is more difficult, changes in the threshold parameter in one 

25 embodiment can be observed when changes occur in either the shape or scale parameter 
without a corresponding change in the other. As will be discussed below with respect to 
Figure 7, control charts may be created in either a logarithmic scale or in an "antilog" 
scale, which is the same scale in which the original moisture content data was taken. 
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5 Figure 6 illustrates a block diagram which describes, in one embodiment, the 

process of calculating estimates for fi and cr. This process can be perforaied as the 
process of blocks 512 and 514 of Figure 5. Ristea and Maness, 2001, demonstrates that 
this process is equivalent to performing the calculations given above in (3) and (4). The 
illustrated embodiment assumes that the known good data set comprises m samples with 
10 each sample comprising n moisture content measurements. In a preferred embodiment, m 
is at least 20, while /j is 5 or greater. At block 602, each sample x^,x^ ,K ,x„ is 
transfomied to normality by first subtracting the estimated threshold parameter and then 
taking the natural logarithm, giving a nomialized sample 7, , 3^2 , by the equation: 

\n{x,-e)^y, (5) 
1 5 Next, at block 604, calculations are made for each of the m samples. In block 

606, for each sample a mean normaUzed value, y , is calculated by taking the mean of 
each of the normahzed measurements: 

-^A±2VtK±A (6) 
n 

Similarly, in block 608, the standard deviation, S, of each sample is calculated. 
20 Next, at block 610, the parameters /u and crare calculated from the y 's and S's by taking 
the mean over every sample. Thus, if 3?,. and are the mean and standard deviation of 
the sample, respectively, and there are m samples, then: 

ji^ = (yi+>^2+K:+3;„)//« (7) 
5=(5,+5'2+K+5,)/w (8) 

25 Because 3; is an estimate of the population mean of the Normal variable, it is also 

an estimate of the scale parameter //. Additionally, because of the way y is created by 
mathematical manipulation of the logarithmically-normalized moisture content 
• measurements, >/ is in fact the geometric mean of those measurements, after taking 
threshold scaling into account. As mentioned above, the geometric mean is known to be 
30 a better indicator of the central tendency of Lognoraial distributed variables. 
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5 S , however, is not an unbiased estimator of the Lognormal shape parameter. In 

order to provide an unbiased estimate of a; at block 616, in the preferred embodiment, 
S is divided by a constant, , the value oif which depends on the sample size. The values 
of C4 for different sized samples can be found in the Hterature. Thus, an estimate for a is 
provided by 5 / C4 in block 616 and the process of estimating the Lognormal parameters 
10 is completed. 

Figure 7 illustrates a block diagram of, in one embodiment, a process used to 
create control charts for the scale and shape Lognormal parameters. Block 516 represents 
the process of developing the charts themselves before new data is plotted on them. 
Beginning with block 704, a sample size n is established for the all new data that will be 

1 5 measured as part of the SPC process. The sample size must be estabUshed because the 
control limits for the scale parameter are dependent on the sample size, as will be shown 
below. Next, at block 706, the center lines of the control charts are chosen. The center 
Imes are sunply /i and<T , as the control charts will be measuring changes in scale and 
shape of the Lognormal distribution. Because the estimated fianda parameters have 

20 been calculated in the process illustrated in Figure 6, the center Imes are ahready known 
to be y and 5 / C4 for the scale and shape control charts^ respectively. Next, in block 708, 
the control limits for each of the control charts are calculated. 

For the scale chart, control limits are established is the traditional fashion for a 
control chart of averages. Thus, in block 710, in one embodiment, conventional "3- 

25 sigma" Ihnits are established around the center line y . The plotted data points will be 
averagesof the transformed measurement samples, i.e. y. The sample means y is 

normally distributed and has standard deviation &/y/n. Hence the Upper Center Line, 
Center. Line, and Lower Center Line of the control chart are: 
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UCL =3; + 3* 




Center Line = y 



(9) 



S 




For the shape chart, because the sample standard deviation is not normally 
distributed, 3-sigma limits cannot be used. However, it is known that the variance of a 
Normal variable such as 7 is chi-square distributed. Thus, m one embodiment, at block 
712 probability limits are used to calculate the control limits. In one embodhnent, this is 
done by the following calculations: 



While the calculations given in (9) and (1 0) create control charts in one 
embodiment, because of the logarithmic transformation done on the data, the scale of the 
charts themselves has no units. Thus, it may not provide the detail of information desired 
when monitoring the wood drying process. Thus, in one embodiment, at block 714, the 
limits and center line are scaled to be put back into the original scale of measurement. In 
one embodiment this is done by taking the antilog of the limits or center line through 
exponentiation and then, in the case of the scale parameter, adding back the estimated 
threshold. The threshold does not affect the standard deviation and so does not need to be 
added back for the shape parameter chart. Thus the calculations in (9) become, in scaled 
form: 




Center Line - — . 

^4 



(10) 
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LCN = e + GKp 
and the calculations of (10) become: 



UCL = 0 + exAy + 3* ^ 



Center Line = 9 + exp^) 
y-3* 



(11) 



C/CL = exp 
Center Line = exp 

iCiV=.exp 



5" \X^al2,n-\ ^ 
«-l J 



^5 



(12) 



{cj n-l J 



It will be demonstrated later that changing the scale in this manner does not 
qualitatively affect the general trends of the control charts, and thus the procedure of 

10 block 714 can be performed or not performed depending on the need of whoever or 
whatever is monitoring the control charts. Finally, at block 716, new data is plotted on 
the chart in order to monitor the drying process. 

The process of block 716 is described in further detail with respect to Figure 8. 
Figure 8 shows, in one embodiment, the process of taking new data, transforming it, and 

1 5 plotting it on control charts so as to monitor the drying procedure. Startmg at block 802, 
a new set of wood moisture content data is gathered. In a preferred embodiment, the 
gathering is done randomly and independently. In another preferred embodiment, the 
data is gathered from the same type of wood dried under the same conditions as that used 
to make the control charts. In yet another preferred embodiment, the new data is 

20 gathered in samples of the same size as was established at block 704. This is done to 

ensure that the proper calculations were performed for the control lines against which the . 
new data will be measured. Continuing to block 808, the data is transformed so that it 
may be plotted on the chart. In block 810, each measurement is transformed from the 
original moisture content into a normalized logarithmic form, as was done for the known 
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5 good set of data in block 602. Next, the data is manipulated by sample so that it may be 
plotted. In block 814, the mean of each transfqraied sample in one embodiment is 
calculated if the chart being plotted on is a scale chart.^ This is done by taking the mean 
of all of the transforaied values for a single sample, creating a y value that can be plotted, 
similarly to the process described above in block 606. Thus, for each 
1 0 sample ;c, , ,K , of size w, each plotted value y is computed by: 

y = i*Xhi(x,-.^) (13) 

In block 816, in another embodiment the standard deviation of each sample of 
transformed data is calculated to detemiine an S value which can be plotted on a shape 
chart. Similarly to the discussion above with respect to Figure 7, the data determined 

15 now is unitless; if the same scale as the original moisture content data is desired, the 

calculated scaled antilog of the J? and S values mixst be found. Thus, in one embodiment, 
the y and S values are changed to the original scale. This is done, similarly to above, by 
exponentiating the values, and, if the chart is a scale chart, adding the threshold 
parameter back in. For example, for the scale chart the value to be plotted can be 

20 determined by: 

y«miog =^ + exp(J?,,g) 

Because y is the logarithmically-normalized mean of the original sample, y is 
related to the geometric mean of the original sample measurements. Thus, in an alternate 
embodiment, the original-scale plot values can be found for a sample by: 
25 y^^^=0■\rgeomean{x,-9) ^^^^ 

Where geomean is the geometric mean. Thus (15) can be peformed directly from 
the measurements by: 

i 

^ ' . ^ (16) 
Next, at block 820, the value is ready to be plotted on its control chart. Then, at 
30 block 822 in one embodiment, it is detemiined if the kita is in ah out-of-control state. In 
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5 one embodiment detennination that the kibi is out-of-control is done because a single 
J; or iS value falls out side of the set control limits. In another, a minimum set number of 
data points are required to fall outside the control limits before an out-of-control state is 
assumed. In another, an out-of-control state may be inferred by the proximity of data 
values to the control limits, regardless of whether and exceed the actual limits. 

10 Figures 9a-d demonstrate examples of control charts for scale and shape, and 

illustrate the similarities between unitless charts and charts in the original moisture 
content units. In one embodiment these charts are displayed as new data is transformed 
and compared to the control limits and center lines computed above. In another 
embodiment the charts are not displayed, and other indications are made about the 

IS relation between the transformed data and the control limits and center lines. The charts 
of Figures 9a and 9b represent the same data for the scale parameter, with the chart of 
Figure 9b having been transformed back into original unit scale. The charts of Figures 9a 
and 9b are known as "y bar" and "antilog y bar" respectively. There is a symmetry shift 
between the two charts; the chart of Figure 9a is symmetrical around its center line due to . 

20 the calculations in (9) while the chart in Figure 9b is asymmetrical due to the 

exponentiation. However a graphical comparison between the two charts illustrates that 
trends from sample to sample are similarly represented in both the unitless and original 
scale charts. Similar tendencies can be observed in the unitless and original scale charts 
for the shape parameter in Figures 9c and 9d. Thus, flie act of scaling discussed above 

25 with respect to Figures 7 and 8 is reasonable and in various embodiments can provide 
similar monitoring abilities as the unitless charts. 

The method described above provides an improved process for statistically 
analyzing wood moisture content data in order to monitor drying kilns. The method 
produces control charts which can be easily monitored, by computer or by human, in 

30 order to deteraiine if a given drying process is in-control. In view of the many possible 
embodiments to which the principles of our invention may be applied, we claim as our 
invention all such embodiments as may come within the scope and spirit of the following 
claims and equivalents thereto. 



